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Abstract 

We consider planar lattice walks that start from a prescribed position, take their 
steps in a given finite subset of and always stay in the quadrant x > 0,y > 0. We 
first give a criterion which guarantees that the length generating function of these walks 
is D-finite, that is, satisfies a linear differential equation with polynomial coefficients. 
This criterion applies, among others, to the ordinary square lattice walks. Then, we 
prove that walks that start from (1,1), take their steps in {(2, —1), (—1, 2)} and stay 
in the first quadrant have a non-D-finite generating function. Our proof relies on a 
functional equation satisfied by this generating function, and on elementary complex 
analysis. 


1 Introduction 


The enumeration of lattice walks is one of the most venerable topics in enumerative com¬ 
binatorics, which has numerous applications in probabilities [^, |^. These walks take 


their steps in a hnite subset © of and might be constrained in various ways. One can 
only cite a small percentage of the relevant litterature, which dates back at least to the 
next-to-last century III, O O 120, |20, BDI . Many recent publications show that the topic is 
still active i, 0,1, 0 as 11^ 

After the solution of many explicit problems, certain patterns have emerged, and a more 
recent trend consists in developing methods that are valid for generic sets of steps. Special 
attention is being paid to the nature of the generating function of the walks under con¬ 
sideration. For instance, the generating function for unconstrained walks on the line Z is 
rational, while the generating function for walks constrained to stay in the half-line N is 
always algebraic [Q]. This result has often been described in terms of partially directed 2- 
dimensional walks conhned in a quadrant (or generalized Dyck walks [0, 0, 0 0, HI), 
but is, essentially, of a 1-dimensional nature. 


^Partially supported by the INRIA, via the cooperative research action Alcophys. 
^Partially supported by MZT RS under grant PO-0511-0101. 
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Similar questions can be addressed for real 2-dimensional walks. Again, the generating 
function for unconstrained walks starting from a given point is clearly rational. Moreover, 
the argument used for 1-dimensional walks conhned in M can be recycled to prove that the 
generating function for the walks that stay in the half-plane a; > 0 is always algebraic. What 
about doubly-restricted walks, that is, walks that are conhned in the quadrant x > 0,y > 07 
It would be satisfactory if the hierarchy unconstrained walks / mono-constrained walks / 
bi-constrained walks could match the classical hierarchy of generating functions: rational 
series, algebraic series, D-£nite seriesQ (also called holonomic series). A rapid inspection of 
the most standard cases only corroborates this hope. For instance, the generating function 
for walks on the square lattice (with North, East, South and West steps) that start from the 
origin and stay in the hrst quadrant is 

^ J\[m/2\J{[n/2\) ^^{ln/2\J {\n/2]J ’ 

which is a D-£nite series. The hrst expression comes from the fact that these walks are 
shuffles of two prehxes of Dyck walks, and the Chu-Vandermonde identity transforms it into 
the second simpler expression. 



Figure 1: A walk on the diagonal square lattice conhned in the hrst quadrant. 


The case of the diagonal square lattice, where the steps are North-East, South-East, 
North-West and South-West (Figure |I[) is even simpler: by projecting the walks on the x- 
and y-a.xes, we obtain two decoupled prehxes of Dyck paths, so that the generating function 
for walks in the hrst quadrant is now 


E 

n>0 


n 

[«/2J 


f 


another D-hnite series. In both cases, the number of n-step walks can be shown to grow 
asymptotically like 4’^/n, which prevents the corresponding generating function from being 
algebraic (see for the possible asymptotic behaviours of coefflcients of algebraic series). 


series F{t) is D-finite if it satisfies a linear differential equation with polynomial coefficients in t. 
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In Section of this paper, we shall generalize this result by proving that, if the set of steps 
© is symmetric with respect to the x-axis and satishes a small height variation condition, 
then the generating function for walks with steps in ©, starting from any given point {ioijo), 
is D-£nite. This result covers the above two cases. 

However, and most importantly, we shall also prove in Section that this holonomy 
result does not hold for any set of steps: walks that start from (1,1), take their steps in 
© = {(2, —1), (—1, 2)} and always stay in the hrst quadrant have a non-D-£nite generating 
function. The central point of our proof is the study of a series G'(x), dehned by an equation 
of the form 

G(x) + G(^(x)) = H(x), 


where ^(x) and A{x) are explicit algebraic series in x. We consider the solution G(x) to 
this equation as a function of a complex variable x, and prove that it has infinitely many 
singularities, which prevents it from being D-£nite. Hence our proof is based on complex 
analysis. To our knowledge, there is no classihcation of the solutions to this type of equation. 
In some very specihc cases (like ^(x) = x^ or ,^(x) = cx) some hypertranscendence results^ 
have been obtained, either by some ad hoc methods [|T^ |^, ^ , or via general results about 
“very” lacunary series . 

These two sections raise the question of a classihcation of the sets © according to the 
nature of the generating function for walks in a quadrant that take their steps in ©. Let 
us mention that some sets of steps, like © = {(1,1), (0, —1), (—1,0)} yield, for non-trivial 
reasons, algebraic generating functions 

Finally, in Section we say a few words about the closely related topic of multidimen¬ 
sional linear recurrences with constant coefficients, and prove a non-holonomy result that 
was announced (but not proven) in H. 


Let us conclude this introduction with a few more formal dehnitions on walks and power 
series. 

Let © be a hnite subset of A walk with steps in © is a hnite sequence w = 
{wo,Wi,... ,Wn) of vertices of if such that Wi — G © for 1 < i < n. The number 
of steps, n, is the length of w. The starting point of w is wq, and its endpoint is Wn- The 
complete generating function for a set A of walks starting from a given point wo is the series 

n>0 i,j£Z 

where aij{n) is the number of walks of A that have length n and end at (i, j). This series is 
a formal power series in t whose coefficients are polynomials in x,yAlxA/u- We shall often 
denote x = 1/x and y = 1/y. The length generating function for walks of A is simply 

A{t) = a{n)A 

n>0 


where a{n) is the number of walks of A that have length n. Note that A{f) = A(l, l;t). 

series F{x) is hypertranscendental if it does not satisfy any polynomial differential equation of the 
type P{x, F{x), F'{x),..., F^'^'>{x)) = 0, where P is a polynomial. 


3 






Given a ring L and k indeterminates Xi,...,Xk, we denote by L[xi,..., 
(resp. L[[xi,... ,Xfc]]) the ring of polynomials (resp. formal power series) in xi,... ,Xk with 
coefficients in L. If L is a held, we denote by L(xi,..., x^) the held of rational functions in 
xi,... ,Xk with coefficients in L. 

Assume L is a held. A series F in L[[xi,... ,Xfc]] is rational if there exist polynomials 
P and Q in L[xi,... ,Xfc], with Q ^ 0, such that QF = P. It is algebraic (over the held 
L(xi,..., Xfc)) if there exists a non-trivial polynomial P with coefficients in L such that 
P{F, xi,..., Xfc) = 0. The sum and product of algebraic series is algebraic. 

The series F is D-finite (or holonomic) if the partial derivatives of F span a hnite di¬ 
mensional vector space over the held L(xi,... ,Xfc) (this vector space is a subspace of the 
fraction held of L[[xi,... ,Xfc]]); see for the one-variable case, and otherwise. In 

other words, for 1 < i < A:, the series F satishes a non-trivial partial diherential equation of 
the form 


E = 0. 


where P^^i is a polynomial in the Xj. Any algebraic series is holonomic. The sum and product 
of two holonomic series is still holonomic. The specializations of a holonomic series (obtained 
by giving values from L to some of the variables) are holonomic, if well-dehned. Moreover, if 
F is an algebraic series and G{t) is a holonomic series of one variable, then the substitution 
G{F) (if well-dehned) is holonomic ||25|, Prop. 2.3]. 


2 A sufficient condition for holonomy 

Let © be a hnite subset of We say that & is symmetric with respect to the x-axis if 

(i,j) G 6 ^ (b-j) e 6. 

We say that © has small height variation if 

(i, j) G © ^ IjI < 1. 

The usual square lattice steps satisfy these two conditions. So do the steps of the diagonal 
square lattice (Figure IH). 

Theorem 1 Let & be a finite subset of 1? that is symmetric with respect to the x-axis and 
has small height variations. Let {io,jo) ^ Then the length generating function for walks 
that start from {io,jo), take their steps in © and stay in the first guadrant is D-finite. 

We shall need the following preliminary result, which does not require any property on ©. 

Proposition 2 Let & be a finite subset oflfi. Let {io,jo) G N x Z. Then the complete 
generating function for walks that start from (zq, jo)? tO'ke their steps in © and stay in the 
right half-plane x >0 is algebraic. 
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Proof. This result is basically of a 1-dimensional nature: projecting the walks on the x-axis 
reduces it to the enumeration of walks on the half-line N, starting from in which each 
step of size i is weighted by a Laurent polynomial in y. 

Y. y'- 

The weight of a walk is taken to be the product of the weights of its steps. The material 
can be readily extended to weighted paths, where the weights belong to any 
held of characteristic 0 (here, Q{y)), and we conclude that the complete generating function 
for walks in the right-half plane is algebraic over Q{x,y,t). 


of 


or 


h6 


Note. The above proposition is very close to the results on partially directed 2-dimensional 
walks mentioned in the introduction. However, in most references, the authors only consider 
the case where all steps are directed; that is, i > 0 for all step {i,j) G ©. 

Proof of Theorem |l]. Let Q denote the set of walks that start from (io,jo), take their steps 
in © and stay in the hrst quadrant. We shall prove that walks of Q are, roughly speaking, 
“equivalent to” walks in the right half-plane ending on the x-axis. 

We claim that it suffices to prove that the subset of Q consisting of the walks that hit the 
x-axis at some point has a D-hnite generating function. Indeed, the remaining walks are, by 
a vertical translation, in one-to-one correspondence with walks that start from (io,jo — k), 
for some k G [1, jo]) stay in the hrst quadrant and hit the x-axis. 

Let us hrst focus on the set of walks in the hrst quadrant that start from (io; jo)) hit 
the x-axis and end at an even ordinate. These walks are in bijection with the set TLq of walks 
that start from (io) jo)) stay in the right half-plane and end on the x-axis. This bijection, 
illustrated by Figure |^, is a mere adaptation of a classical 1-dimensional correspondence, 
which establishes that Dyck prehxes ending at an even ordinate are equivalent to bilateral 
Dyck walks (see the Catalan factorisation in |^.) Starting from a walk w of Q^\ ending 
at level 2k, we denote by Si,..., the steps that follow the last visit of tn to a point of level 
(ordinate) 0,..., /c — 1. Replacing these k steps by their symmetric steps with respect to the 
x-axis yields a walk w that belongs to Ho- The ordinate of the lowest point(s) visited by w 
is —k. Conversely, the steps si,... ,Sk of w that we have to flip back to recover w are the 
hrst steps of w that lead to level — j, for l<J<k. 

Let H{x, y; t) be the complete generating function for walks in the right half-plane, start¬ 
ing from (io) jo)- By Proposition H, this series is algebraic, hence D-£nite. The length 
generating function for walks of Tjo is obtained by extracting the coefficient of y^ in the 
generating function H{l,y;t). Extracting the constant term of a D-hnite series is known to 
give another D-hnite series [^: thus the generating function for walks of Tio, and hence of 
Q!'e\ is D-hnite. 

A similar argument holds for the set of walks in the hrst quadrant that start from 
(zq, jo)) hit the x-axis and end at an odd ordinate: they are in one-to-one correspondence 
with the set 7j_i of walks that start from {io,jo), stay in the right half-plane and end at 
level —1. The generating function for walks of 7j_i is the coefficient of y^ in yH{l, y; t), and 
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Figure 2: The bijection between some walks in the quadrant and walks in the right half-plane 
ending on the x-axis. 

hence is D-finite. Given that the sum of D-finite series is D-finite, this concludes the proof 
of Theorem ^ 

■ 

Remark. A slightly stronger result is proved in via a functional equation approach: 
under the assumptions of Theorem the complete generating function Q{x,y]t) for walks 
in the quadrant, counted by their length and coordinates of the endpoint, is also D-£nite. 

3 The knight walk is not holonomic 


3.1 The main result 

We study walks that start at (1,1), take their steps in {(—1, 2), (2, —1)} and stay in the first 
quadrant. We call them knight walks, since their steps correspond to two of the knight moves 
on a chessboard (Figure |). We note that a walk ending at {i,j) always has i + j — 2 steps: 
hence the information contained in the complete generating function is actually already 
contained in the following bivariate series: 

Q{x,y)= 

*>0j>0 

where Qij denotes the number of knight walks ending at 

The coefficients Qij satisfy the following recurrence relation: 

r 0 if i < 0 or j < 0, 

Qi,j = w if * = i = 1, (1) 

[ Qi+ij-2 + Qi-2,j+i otherwise. 

This recurrence allows us to compute the numbers Qij inductively, for instance diagonal by 
diagonal. The first few values are given in Table |l| below, in which the zero entries are left 
out. The non-zero entries he on the lines i = j mod 3. 
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Figure 3: A knight walk. 


Applying the transformation (i,j) —> ((2z + j)/3, {i + 2j)/3) shows that Qij is also the 
number of walks made of North and East steps, that start from (1,1), end at {{2i + j)/3, {i + 
2j)/3) and always stay above the line 2y = x and below the line 2x = y (see Figure @). 
One thus recognizes a problem already studied by Gessel |T^, Section 4], But the results 
presented here are stronger than those of 


ISfl- Ignoring the two boundary lines gives the 


following simple bound: 


In particular, 


Qi,j — 


i+j -2 
(2z+j-3)/3 


Q3ifl ^ 


3i-2 
i - 1 


This bound could be sharpened by counting walks that stay above the line 2y = x 
e.g. [|^ or 1^, Example 4]), but the above bound will be enough for our purpose. 


( 2 ) 


see 


Figure 4: Another description of the knight walk. 
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j 

T 

12 

11 

10 

9 


24 

24 

6 


7 

6 

5 

4 

3 

2 

1 

0 


2 

1 


0 


6 

2 

1 

I 


108 

84 

24 

18 

6 

4 

2 

1 

"I ^ 


312 

60 

12 

2 

T” 


204 

36 

6 


5 


720 

408 

120 

60 

18 

6 

2 

“6 f 


1440 


204 

24 


720 

312 

84 

24 

6 

9 10 


108 

24 

Tl 12 


Table 1; The number Qij of knight walks ending at 


By summing the recurrence relation we obtain: 

y) = ^ (Qi+i,i-2 + Qi-2,j+i)x'-y^ 

= xy + y^/x{Q{x, y) - Q(0, y)) + x'^/y{Q{x, y) - Q{x, 0)), 

that is 

{xy -x^ - y^)Q{x, y) = x^y'^ - G{x) - G{y) (3) 

where 

G{x) = x^'^Qifix'= x^Q{x,Gi) (4) 

i>0 

counts knight walks ending on the x-axis. (We have used the symmetry of the problem in 
X and y.) Note that the length generating function for all knight walks is t~‘^Q{t,t). The 
above equation, combined with the elementary properties of D-finite series, imply that the 
following three statements are equivalent: 

— the bivariate generating function Q{x, y) for knight walks is D-finite, 

— the generating function Q{x, 0) for knight walks ending on the x-axis is D-finite, 

— the length generating function t~‘^Q{t,t) for knight walks is D-finite. 

Our main result asserts that none of these statements hold. 

Theorem 3 The length generating function for walks that start from (1,1), take their steps 
in {(—1, 2), (2, —1)} and always stay in the first guadrant (knight walks) is not D-finite. Nor 
is the generating function for knight walks that end on the x-axis. 

Looking at Table |l|, one might still have some hope that the numbers Qtj are not so bad. In 
particular, they seem to have small prime factors. This pattern actually does not go on, and, 
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in case there would still be a doubt, the following proposition reinforces the non-holonomic 
character of these numbers. 

Proposition 4 The length generating function for knight walks ending on the main diagonal 
X = y is not D-finite. 

3.2 The kernel method 

The so-called kernel method solves completely the functional equation (^. This method has 
been around since, at least, the 70’s, and is currently the subject of a certain rebirth (see 
the references in iSi)- 

Applied to our equation, this method consists in coupling the variables x and y so as to 
cancel the kernel xy — x^ — y^] this yields the missing information about the series G{x). 
More precisely, let be the unique formal power series in x satisfying 

xf — x^ — = 0. 

The Lagrange inversion formula, applied to f{x)/x, provides an explicit expression for ^(x): 



Replacing y by in gives a functional equation that dehnes the power series G{x)-. 

G{x) + G{f{x))=x^a^f. ( 6 ) 

Indeed we obtain, after iterating this equation inhnitely many times: 

G(x) = 5^(-l)' (7) 

i>0 

where ^ o • ■ • o ^ is the Ah iterate of f. Note that i^''\x) = so that the sum is 

convergent in the ring C[[a:]]. Replacing G{x) by the above explicit value in would give 
an expression for Q{x,y). 

However, we shall not exploit these expressions, but rather the functional equation (j^), 
to prove that G{x), hence Q{x,y), is not D-£nite. Our proof will be of an analytic nature. 
The idea is to consider G{x) as a function of a complex variable x and study its singularities. 
Using functional equations like (|^), we shall build new singularities of G from old ones - and 
end up with inhnitely many singularities, thus proving that G cannot be holonomic. 

However, even though (|) dehnes G{x) uniquely, this equation itself is not sufficient for 
our purpose: we shall introduce the other two roots and ^ 2 ( 2 ^) of the kernel, and the 
corresponding analogues of (|), to obtain enough singularities. 


9 



3.3 The roots of = xy 

As a polynomial in y, the kernel xy — x^ — y^ has three roots. Only one of them, given by (|^), 
is a power series in x. We shall denote it, from now on, by 

.^o(a;) = + 3x® + + 55x^^ + 273x^’^ + ■ ■ ■ 


The other two roots are power series in yT, and their expansion can be compnted indnctively: 


^i(^) = 

= 


X 


-X vx 


X 

Y 


X 


105 

105 


x^\/x 


Q rp^ 

J-' ^3 r~ -LKJU 


Of conrse, ^ 2 (x) is derived from .^i(x) by replacing ^/x by —^/x. Gnided by the above 
expressions, let ns write 


,^i(x) = +^/x■^|J{x) - (j){x) (8) 

^2{x) = -\/x'?/>(x) - 0(x) 

where "0 and 0 are formal power series in x. As the three roots snm to zero, one has 


0(x) = ^o(x)/2. 


In order to compnte the coefficients of '0, we shall nse again the Lagrange inversion formnla 
(LIF). Let x(x) be dehned by 


6Y) 


\fx 

+ X{x) 


Then 

X = x^/x{l + xf^‘^ , 


so that the LIF gives 



Separating the even/odd valnes of n gives explicit expressions of '0(x) and 0(x). The following 
lemma snmmarizes the resnlts thns obtained. 


Lemma 5 Let ^ and -0 be the following power series in x: 




x‘^Y. 

m>0 


/3m\ 
2m + 1 \ m / 


i-E 

m>l 


m!(6m)! x^™' 

(6m — l)(2m)!2(3m)! IG’" 
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Then the three roots of = xy are 


fo{.x) = f 

^i(a;) = ^/x'^l)-f/2 

f 2 {x) = -y/xf)-f/2. 

Both f{x) and V’(a^) have radius of convergence Xc = 4^/^/3. 

The last statement is obtained using the Stirling formula. Note that one also has the following 
closed form expressions: 

^ fx . { arcsin ((3a:)3/V2)^ 

{(U = -3- )• 

arcsin ((3x 
3 




We shall now consider the ^i(x) as functions of a complex variable x. We choose a 
determination of the square root that coincides with the usual determination on R’*': 

V re*® = y/re^^^"^ for — tt < 6 < tt. 

Figure ^ shows the (real) values of the functions fi for real values of x. When x is positive, 
the plots show, from bottom to top, ^ 2 {x), ^o{x) and ^i(x). When x is negative, the only- 
real branch is ^o- 



Figure 5: The real values of fi{x) for a real x. 

The implicit function theorem implies that the singularities of the functions fi are to be 
found among the complex numbers x such that the pair {x,y) = {x,f,i{x)) satisfies 3y^ = x 
(and, of course, x^ + y^ = xy)] that is, 

{x,y) G {(0,0), {xc,yc), (jxcfyc), ifxcjyc)} (9) 

where Xc = 4^/^/3, yc = 2^/^/3 and j = exp(2i7r/3). A more detailed investigation gives the 
following result. 
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Lemma 6 The singularities of the functions fi are given by: 


Sing{fo) = {xcjxcfxc}, Sing{fi) = {0,Xc} and Sing{f 2 ) = {Ojxcfxc}. 


In particular, ^2 is not singular at Xc, as suggested by Figure ||. All the above singularities 
are of the sguare root type. 


Proof. For each of the values of x given by (^, we hrst have to compute the values ii{x), 
0 < i < 2, in order to determine which pairs {x,fi{x)) are actually critical, and then, to 
check the existence and nature of the singularity. 

At a; = 0, all the ^i{x) are zero. The explicit expansion of Lemma |^, combined with 
the fact that ip and f have a positive radius of convergence, shows that only ,^1 and ^,2 are 
singular — their singularity being obviously of a square root type. 

When X = x^, factoring the polynomial x^ + y^ — xy shows that the multiset 
{■^ 0 ( 3 ^ 0 );'Ci(3^c),'C 2 (a^c)} equals {|/c, 2 /c 5 —2j/c}. The functions fi{x) are real and continuous 
on [0, xf\, and the expansions of Lemma | show that ^q^x) and are positive as x —>• 0+, 
while ^ 2 ( 2 ^) is negative. As the fi can only vanish at 0, this sign pattern must go on until Xc, 
so that 

io{xc) = = Vc and i 2 {xc) = -2yc. (10) 

In view of (H), Xc cannot be a singularity of ^ 2 - Now a local expansion of x^ + y^ — xy around 
{xcVc) shows that as x approaches xf, 


^ 0,1 (a^) 



a/^ 


x(l + 0(1)), 


( 11 ) 


which conhrms that the singularity of .^0 and .^1 at Xc is of the square root type. 
The values (|l0|) and Lemma |] imply in particular that 


^c'pi^Xc) = 

Using this result, the fact that ip and f are essentially functions of x^, and the values 

V2 = -f' and = -T 

we now compute 

^o(ja;c) = ^o(j^a;c) = jVc 

i\{3x<P) = -23^yc = -2fy^ 

^(ja^c) = d^Vc ^2{fxc) = jyc 

so that cannot be singular at jxc or j^Xc. Finally, local expansions of x^ + y^ — xy conhrm 
as above the existence of square roots singularities of ^0 and ^2 at jXc and j^Xc. 


Each of the functions fi has a unique analytic continuation on any simply connected 
domain avoiding Sing{^i), for instance the domain obtained by removing from C the four 
half-lines of Figure 
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Figure 6: A domain on which all the functions are holomorphic. 

Observe that the series and ^2 can also be substituted for y in the functional equa¬ 
tion (U); thus for i G {0,1, 2}, the following equation holds 

G(a;) + G(6(x)) = (12) 

at least as an identity between power series in x or ^/x. 

We end this subsection with a lemma that will be useful to build large singularities from 
small ones. 

Lemma 7 Let x G C, x ^ 0. Then one of the roots of x^ + — xy = 0 has modulus larger 

than |x|. 

Proof. Let yo,yi,y 2 denote the three roots, and assume none of then has modulus larger 
than |a:|. The relation yoyiy 2 = —x^ forces ||/o| = |?/i| = I 1 / 2 I = |a;|. Then, the relation 
2/0 + 2/1 + 2/2 = 0 implies { 1 / 1 , 1/2} = {jyo,fyo}. Finally, the relation yoyi + 2/02/2 + 2/i2/2 = -x 
yields x = 0. ■ 

3.4 The series G{x) is not D-finite 

We now turn our attention to the series G{x) dehned by (!)• 

Proposition 8 The series G{x) has radius of convergence x^ = 4^/^/3. It is singular at Xc, 
and, as X —>■ xf, 

G{x) = A — i?A/W-T7^(l + 0 ( 1 )), 
where A and B are non-zero real numbers. 

Proof. We start from the functional equation dehning G\ 

G(x) + G(^o(a:)) = x^fo(xf- 

Recall that G(x) and GG) have nonnegative coefficients (G(x) counts walks, and .^ 0 ( 2 ^) = 
f(x) is given explicitly by (H)). Hence G(^o( 3 i)) and x^fo(xy also have nonnegative coeffi¬ 
cients. The radius of convergence of ^0(2;) being Xc, the series G(x) and G(^o(a^)) have radius 
at least Xc- The fact that G{x) has radius at least Xc can also be directly derived from the 
upper bound (0). 
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For |x| < Xc, one has |^o(2^)| < ^o(3^c) = He, and Uc is smaller than Xc- This implies that 
the above functional equation also holds as an identity between analytic functions of x in the 
disk |x| < Xc- As X approaches Xc inside this disk, a local expansion gives, thanks to (pd^, 

G(x) = xlvl - G{yc) - - x (2x^2/c - G'{yc)) + 0{xc - x). 

The upper bound (|^), combined with (^), yields 

i>l ' 

This sum can be evaluated numerically, and is found to be smaller than 0.16. Comparing 
with 2x^|/c ~ 0.23 gives the announced result. 


In view of Proposition ||, Eq. (|^ now holds as an identity between functions of x, as 
long as |x| < Xc and |Cj(x)| < Xc (and x ^ if i = 1 or 2). 

Proposition 9 The series G{x) that counts knight walks ending on the x-axis is not D-finite. 


Proof. Assume G is D-£nite. Then, it has a finite number of singularities, and has a unique 
analytic continuation on any simply connected domain avoiding these singularities. For 
i G {0,1,2}, the series G{fi) are also D-£nite (since fi is algebraic). Moreover, by analytic 
continuation of (|I^), 


G(x) + G(ii(x)) = x^(,(x)- 


(13) 


for all X in any domain where all our functions are analytically defined. Let us consider the 
identity (|^) for i = 2. As x approaches Xc, the function G(x) becomes singular, while 
does not: this shows that ^ 2 {xc) = —2.yc is a singularity of G. 

Now let Xg be a singularity of G of maximal modulus. According to Lemma ^ there 
exists i such that |^j(xs)| > |xs|. By assumption, G(x) is singular at Xg- But |xs| > 2yc > Xc, 
so that Xg is not a singularity of fi- Hence ([T^) implies that G is singular at fi{xs). But 
|.^i(xs)| > |xs|, which contradicts the maximality of |xs|. 


Remarks 

1. The only property of D-£nite series we have really used is the fact that these series have 
hnitely many singularities. Hence, what we have actually proved is that G cannot have 
finitely many singularities. 

2. The principle of the proof can be applied to other functional equations of the same type. 
We apply it in the next section to a (minor) variation of the knight walks. 

3. As mentioned at the beginning of this section, the same walk problem was studied by 
Gessel in [|T^. This reference contains a formula analogous to (^, and the statement that 


G{x) has radius of convergence 4^/^/3. 


There remains to prove Proposition 
on the diagonal is not D-£nite. 


the generating function for knight walks ending 
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Proof of Proposition |^. The length generating function for walks ending on the diagonal 
is the coefficient of in the series t~‘^Q{tx, tx), where x = 1/x. This series is a power series 
in t whose coefficients are Laurent polynomials in x. Let us denote u = x^ and u = 1/u. 
Then, by Eq. (D, 


t'^Qitx, tx) 


- Sjt^u) - Sjt^u) 
1 — t(u + u) 


(14) 


where the series S{z) is dehned by S^z^) = G{z). Let us convert 1/(1 — t{u + u)) in partial 
fractions of u. We obtain 


1 _ 1 ( 1 ^ 1 _ ^ 

1 — t{u + u) V 1 ~ (t) 1 — uU (t) J 


where 


[/(t) 


1 - Vl - 4t2 
2t 


We can now extract the constant term in x (or in u) from ([T^) : 


[x^]t‘^Q{tx, tx) 


1 

Vl - 4f2 
1 

Vl - 4f2 




S{t^u) 

1 - uU{t) 


{t^ - 2S{t^U{t))) . 



S{t^u) \ 
1 - uU{t)) 


Assume this series is D-£nite: then the series D{t) dehned by D(t) = S(t^U(t)) is D-hnite 
in t too. Let T{s) be the unique power series in s such that T(0) = 0, T'(0) = 1 and 
2^6 _ ^ 4^12 gS _ g fg^gi; D-hnite series are stable by any algebraic substitution 

tells us that D{T{s)) = S{T^U{T)) is D-hnite. But T^U{T) = s^, so that ^(s^) itself, and 
hence S'(s) and G(s), are D-hnite too, which we have proved to be false. 


4 A link with multidimensional linear recurrences with 
constant coefficients 

In p, we considered d-dimensional sequences of complex numbers, denoted a„ = ani,n 2 ,...na: 
dehned by recurrence relations of the following form: 

a„ = ^ Chttn+h for n>s, (15) 

h&H 

where H = {hi, h. 2 ,..., h^} C is the set of shifts, are given nonzero constants, 

and s G is the starting point satisfying s + H ^ N'^. We think of the hi as having mostly 
(but not necessarily only) negative coordinates, and of the point n as depending on the 
points n + hi,n + h 2 ,... ,n + hk as far as the value of an is concerned. A given function ip 
specihes the initial conditions: 

Un = pin) for n > 0, nffs. (16) 
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The convex hull of the set H is assumed not to intersect the first orthant {i.e., n > 0). This 
condition, as shown in [§, guarantees that the numbers can be computed recursively 
using (p!5[). 

The enumeration of walks in a quadrant hts exactly in this framework, with d = 3. 
Indeed, denoting Qij{n) the number of walks that start from a given point (fo, jo), end at 
{i,j) and have length n, we have 


Qi,j (^) 


1 if (i, j,n) = (io, jo,0) 

0 if i < 0 or j < 0 or n < 0 

Qi-h,j-k{n - 1) otherwise, 

ih,k)ee 


where & is the set of steps. A translation of the indices i and j transforms this recursion 
into one of the above type, with d = 3 and H = {(—h, —k, —1) such that (h, k) G 6}. 

The paper |P] mostly dealt with the algebraic nature of the generating function of the 
solution of such recurrences, and we started a classihcation, based on the apex of the re¬ 
currence, dehned as the componentwise maximum of the points in if U {0}. We proved 
that when the initial conditions have rational generating functions and the apex is 0 , the 
generating function of the solution is rational. Next, when the initial conditions have alge¬ 
braic generating functions and the apex has at most one positive coordinate, the generating 
function of the solution is algebraic. 

When the apex has two positive coordinates, and d = 3, our study of walks in a quadrant 
shows that the solution might or might not be D-£nite. For 2-dimensional sequences, the 


“simplest” example with apex (1,1) was introduced in [Q: For i,j > 0, let 


dij — 


_ j fli+lJ-2 + 


1 


if i,j > 2, 

otherwise. 


(17) 


This recurrence is obviously closely connected to the knight walk. Again, it has a unique 


solution whose terms can be computed inductively, 
following array. 

j 

T 

6 115 7 
5 1 1 3 5 10 14 

4 1 1 3 4 6 10 

3 1 1 2 2 4 5 

2 1 1 2 2 3 3 

111111 1 
0 11111 1 


The hrst few values are given in the 


7 

5 

1 

1 


0 1 2 3 4 5 6 


Dehning the generating function 


i,j>2 
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we obtain, by summing the recurrence relation over i,j > 2, the following functional equa¬ 
tion: 

{xy -x^ - y^)A{x, y) = R{x, y) - F{x) - F{y) (18) 

where 


R{x,y) = xy 


1 +?/ 1 + 

1 — X 1 — y J 


and F{x) = ^ = x^A{x, 0). 

i >2 


We have used the symmetry of the problem in i and j. It was hrst proved in that F{x) 
and A{x,y) are irrational. Then, we claimed in p, p. 74] that they are even not D-£nite, 
but without giving a proof. The tools developed above for the knight walk apply perfectly 
to this problem. 


Proposition 10 The generating function A{x, y) of the bivariate seguence Oij defined by 
the recurrence relation dH is not D-finite. Nor is the series A(x, 0). 

Proof. The argument is very close to the knight one. The kernel method hrst gives 


F{x) + FiU^)) = R{x,U^)), (19) 

for any of the three roots of the kernel. 

The only difference with the knight treatment comes from the fact that R{x, ^i{x)) might 
have — and indeed, has — more singularities than They can be determined exactly, but 
we shall only use the following obvious information 

Sing{R{x, ^j(x))) C Sing{^i) U {1} U {x : — x -|- 1 = 0}. (20) 

In particular, all singularities of R{x,fi{x)) have modulus at most 1.33 (an upper bound for 
the modulus of the largest root of x^ — x -1- 1). 

We start from Eq. (|I^), in the case z = 0. As R{x, has radius Xc, we End again that 

the radius of F{x) is at least Xc- Moreover, comparing the recurrence relations (|^) and ([T7|), 
and the corresponding tables, shows that for i,j > 2, one has aij > Qi_ 2 ,j- 2 - This implies 
that the radius of F(x) is bounded from above by the radius of G(x), which was proved to 
be exactly Xc- Hence F(x) has radius Xc- 

Assume F{x) is D-£nite. As in the proof of Proposition |^, we hrst construct large singu¬ 
larities of F (here, large means larger than 1.33). Let us start from xq = Xc = 4^/^/3 ~ 0.53, 
which is a singularity of F (by Pringsheim theorem, a series with nonnegative coefficients 
is singular at its radius of convergence). As it is not a singularity of R{x,f 2 {x)), Eq. (pffi]) 
implies that xi := ^ 2 (a^c) = —2z/c ~ —0.84 is singular for F. 

But xi is singular for none of the R{x,fi{x)). Moreover, 

{^o(a:i), ^i(xi), ^ 2 ( 3 : 1 )} = {xo, -0.26... ± 1.02...z}. 

Using the same trick as above, we see that X 2 ~ —0.26 — 1.02z is singular for F. 

One more step: Among the fi{x 2 ), one is X 3 ~ 0.92 — 1.02z, which has modulus larger 
than 1.33. As X 2 is not singular for any of the R{x,fi{x)), the ^i(x 2 ) are singularities of F. 
In particular, X 3 is singular for F. 

We conclude as above, by considering the largest singularity Xs of F (in modulus), showing 
that one of the ^i{xs) is singular for F and larger than |xs| in modulus. 
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